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SUMMARY 


A method is presented which allows the use of nonlinear section 
lift data in the calculation of the spanwise lift distribution of 
unswept wings with flaps or ailerons. This method is based upon lifting- 
line theory and is an extension to the method described in NACA Rep. 865 . 
The mathematical treatment of the discontinuity in absolute angle of 
attack at the end of the flap or aileron involves the use of a correction 
factor which accounts for the inability of a limited trigonometric series 
to represent adequately the spanwise lift distribution. . A treatment of 
the apparent discontinuity in maximum section lift coefficient is also ~ 
described. In order to minimize the computing time and to illustrate 
the procedures involved, simplified computing forms containing detailed 
examples are given for both symmetrical and asymmetrical lift distribu- 
tions. A few comparisons of calculated characteristics with those 
obtained experimentally are also presented. 


INTRODUCTION 


Unswept-wing characteristics calculated by the method of refer- 
ence 1, in which nonlinear section lift data are used, have been found 
to agree much closer with experimental data in the region of maximum 
lift coefficient than those calculated by methods in which linear section 
lift curves are used. It" appears feasible that the similar use of non- 
linear section data would yield improved results for unswept wings with 
flaps or ailerons. The deflection of a partial-span flap or aileron, 
however, causes discontinuities in the spanwise distribution of the 
absolute angle of attack. If such discontinuities exist, an excessively 
large number of spanwise stations must be considered in order to obtain 
a solution by the method of reference 1 or by any other method in which 
the lift distribution is approximated by a trigonometric series. A 
different method of treatment of the discontinuity is therefore desirable. 
Developed herein is a new method of treatment involving the use of a 
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correction factor which accounts for the inability of a limited trigo- 
nometric series to represent adequately the spanwise lift distribution 
of a wing with partial-span flaps or ailerons deflected. This correc- 
tion factor is obtained with the aid of reference 2 and is used in con- 
junction with the system of multipliers of reference 1 to obtain the 
induced angle of attack from the spanwise lift distribution. The sub- 
sequent calculation of the lift distribution by means of successive 
approximations is similar to that of reference 1. 

The mathematical, treatment of the discontinuity in the spanwise 
distribution of absolute angle of attack is only part of the problem 
involved in calculations for wings with flaps or ailerons. The direct 
use of two-dimensional lift data would indicate a discontinuity in the 
spanwise distribution of maximum section lift coefficient. Obviously, 
the flow about the wing sections near the end of a flap or aileron is 
not two-dimensional. A rational method of obtaining three-dimensional 
section data from the two-dimensional data has therefore been devised 
and is presented herein. This method is substantiated by experimental 
pressure distributions. 

In addition to the presentation of the method of makl ng the calcu- 
lations, simplified computing forms are given and their use is illus- 
trated by detailed examples for both symmetrical and asymmetrical dis- 
tributions. A few comparisons of calculated results with experimental 
data are also given. 

This paper and reference 1 are intended to supplement each other. 
The reader is therefore expected to be reasonably familiar with 
reference 1. • 


SYMBOLS 


As used herein, the term "section" designates the characteristic 
of a section in three-dimensional flow. 

A aspect ratio 

A n coefficients of trigonometric series for lift distribution 

Cjw induced drag coefficient 


Cl wing lift coefficient 

C z rolling -moment coefficient 

C ni induced yawing-moment coefficient 
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E edge-velocity factor for symmetrical part of lift distribution 

E 1 edge-velocity factor for antisymmetrical part of lift 

distribution 

F factor used in altering two-dimensional lift curves 

R Reynolds number 

o ~ ~ , . , /Total pressure - Local static pressure\ 

S pressure coefficient * — : ; £ 

\ Dynamic pressure / 

a Q section lift -curve slope per degree 

b span of wing 

c local chord of wing 

c s root chord 

c mean geometric chord (b/A) 

V 

c^ section induced-drag coefficient 

c 1 section lift coefficient 

Ac 1 stall margin 

c 2 ma x maximum section lift coefficient 


c 7 ^ maximum two-dimensional lift coefficient 
^max^ 


Ac 


Z max 


C Z 


* 


c Zi 


c 


z al 


c n 


increment in maximum lift coefficient due to flap deflection 
(two-dimensional data) 

section lift coefficient at either end of flap or aileron 

section lift coefficient for part of lift distribution 
involving no discontinuity in angle of -attack 

section lift coefficient for part of .lift distribution due to 
discontinuity in angle of attack 

section lift coefficient for additional lift distribution 

section normal -force coefficient 
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2V 

s,t 

t 

c 

y 

y* 

a 

a c 

°e 

a i 

a il 

ai 2 

a *o 

a o 

a s 

% 

^mk 


^mk 
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wing-tip helix angle generated by rolling wing 
intervals used in integration 
airfoil thickness -chord ratio 
spanwise coordinate 

spanwise coordinate at either end of flap or aileron 
angle of attack, degrees 

correction for induced angle of attack, degrees 
effective angle of attack, degrees 
correction for effective angle of attack, degrees 
induced angle of attack, degrees 

induced angle of attack for part of lift distribution 
involving no discontinuity in angle of attack, degrees 

induced angle of attack containing discontinuity, degrees 

angle of attack for zero lift, degrees 

angle of attack for two-dimensional lift curves, degrees 

angle of attack of root section, degrees 

uncorrected induced angle of attack, degrees 

multiplier for induced angle of attack for asymmetrical 
distributions 

multiplier for induced angle of attack for antisymmetrical 
distributions 

magnitude of discontinuity, in absolute and induced angles of 
attack, degrees 

angle of' twist, negative if washout-, degrees 
angle of twist at wing tip, degrees 
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e average angle of twist 

e 5 faired angle of twist due to flap deflection 

e g average angle of twist due to flap deflection 

T) ratio of actual two-dimensional lift-curve slope to 

theoretical value of n^/90 

T) m area multiplier for asymmetrical distributions 

Ti ms area multiplier for symmetrical distributions 

9 = cos - ^ 

b 

9 * = cos -1 ^ 

D 

X taper ratio. (| g t C | gg d ) 

multiplier for induced angle of attack for symmetrical 
distributions 

v m interpolation multiplier 

a m moment multiplier for asymmetrical distributions 

a ma moment multiplier for antisymmetrical distributions 

* used as superscript to denote value at end of flap or aileron 


DEVELOPMENT OF METHOD 
Lift Distribution 


The method involving the use of multipliers to obtain the induced 
angle of attack from the spanwise lift distribution was presented in 
reference 1 . In this development, the lift distribution was approximated 
by a finite trigonometric series 


ff ) • e 


a j nut 
A n sin n — 


n=l 


(1) 
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where cos — = ^ .and m = 1, 2, 3»' • • . r - 1. From the values 
r b 777 

of cjc/b for each value of m, the induced angle of attack could then 
he obtained at the points kit/r by the relation 



where k = 1, 2, 3* • • . r - 1 and (3^ denotes the multipliers for 

asymmetrical distributions. The corresponding multipliers, 

and axe used for symmetrical and anti symmetrical distributions, 

respectively. These multipliers are tabulated in reference 1 for r = 20. 
The method can be used directly so long r.s there is no discontinuity in 
the spanwise distribution of absolute angle of attack caused by the 
deflection of a partial-span flap or aileron. 

Lifting-line theory requires that a discontinuity in the distribu- 
tion of absolute angle of attack must be accompanied by an identical 
discontinuity in the distribution of induced angle of attack in order to 
avoid a discontinuity in the spanwise lift distribution. An analytic 
expression for the lift distribution associated with a discontinuity in 
induced angle of attack is presented in reference 2. The complete lift 
distribution can thus be expressed as the sum of two distributions 



(3) 


where c^c^bb is the distribution due' to a unit discontinuity in the 
induced angle of attack and c^c/b is the remainder of the lift dis- 
tribution. These distributions are illustrated in figure 1. Since no 
discontinuity is associated with the distribution c^c/b, the corre- 
sponding induced angle of attack can be obtained by means of the 
multipliers 



. (4) 


By definition, over the flap span 


a 


^k 


= 6 


(5a) 
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and over the unflapped span 


tt i2 k " ° 


(5b) 


The total induced angle of attack is 


a i k " ' a ii, + a i 2 . 


If, however, the multipliers were used with the total distribution 


(6) 


m=l \ /m m=l ' ' “ , m=l 


^ ( c l 2 c ' 

6 b5 y m p mk 


and a^2 were added to both sides of this equation, the result would be 


a c Z c \ - - 

T“ / p mk + a i 2k + XT * 

m=x “ m=i m k fcl 




/=I 2 C\ 

5| — j Pmk 

. 'm 


r-1 


Rearranging this expression gives 




r-1 


a. 





m=l 




Pmk + 6 


a i 




2_ 

h5 


Pmk 


m 


(7) 


A comparison of equation (l), for a wing with a flap or aileron, with 
equation (2), for a wing without flap or aileron, shows the addition of 
a term which is proportional to the magnitude of the discontinuity and 
acts as a correction factor to account for the inability of a limited 
trigonometric series to represent adequately the spanwise lift distribu- 
tion of a wing with partial-span flaps or ailerons deflected. For 
simplicity, equation (7) may be written as 


a ik s v 6 


a r 


( 8 ) 
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where the uncorrected induced angle of attack is expressed as 


= 2_rrY mk 

m=l 


(£l<* <9) 

and the correction factor per unit discontinuity is given by the equation 


a Cfc “1 2 . 




5 _ 5 ~ i— \ 


(10) 


The distribution cj 2 c/b& given in reference 2 may be expressed in 
the form 


C7 c 
l 2 

b6 


ai 2 


90 


, „ r*k s -> 1 - cos(0 + 0*) , Jt0*sin 0 

(cos 8 - cos e ) log. 1 _ cos(9 , y * - ^ 0 — 


(11) 


for = 1, 0° < 0 < 9*-, for ^ =0, 0~ < 0 < l80°j and 0* 

and 6 are in degrees. This distribution is dependent on only the 
spanwise position of the discontinuity and is independent of aspect 
ratio and taper ratio. N The correction factor per unit discontinuity 
Oc^y/6 given by equation (10) is therefore a function of only the span- 

wise position of the discontinuity. The distribution of a,i 2 /&, for 
which equation (11) applies, is illustrated in figure 2(a). Values of 
c j 2 c/bb corresponding to this type of distribution are presented in 

table I for various values of 2y/b' and 2y*/b. These values are 
plotted against 2y/b in figure 3 for even increments in 2y /b and 

against 2y*/b in figure 4 for even increments . in 0 = cos -1 In 

table II are given values of a^/b corresponding to the values of 
cj 2 c/b& of table I.‘ . 

The distribution of illustrated in figure 2(a) would be 

applicable only for a wing with one outboard flap or aileron deflected. 
For a wing with symmetrical inboard flaps, the distribution would be 
obtained as illustrated in figure 2(b) and the values of Oc/& would 
be obtained in like manner. For example, for a wing with flaps 

‘K’ 

extending from = 0.6 to = -0.6, the values of a ^/b to be 

b D 2 * 

used would be obtained by subtracting the values of a c /& for = 0.6 
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from those for = -0.6. Since the resulting distribution would be 

D 

symmetrical, the values of ctc/5 would also be symmetrical. Similarly, 
an antisymmetrical distribution of cl±2/^ the corresponding values 

of oc/b would be obtained as illustrated in figure 2(c). 


It should be noted that two values of 

* 


ctc/S exist at the end of 

the flap, one corresponding to + 0 (flap side of 2y*/b) and one 

2v* ° 

- 0 (unflapped side). Only the values corre- 
are given in table II. s For the unflapped side. 


corresponding to , 
2v* b 

v + 0 

b 


sponding to 


& /2zl. 0 

b 


/a r 


%l + 0 

D 


- 1 


These values have no practical significance unless the discontinuity 

Py kjt 

occurs at one of the stations = cos. — . At these stations, either 

b r * 

of the two values may be used so long as the value is used with the 
proper section lift curve. . 


In the calculation of the induced angle of attack by equation (8), 
the_ values of oc^/6 must be multiplied by the magnitude of the discon- 
tinuity 5. The value of 6 to be used is obtained from the section lift 
curves at the lift coefficient C 2 . If the discontinuity occurs at one 

of the spanwise stations ^ = cos the value of (cje/b)* is obtained as 

one of the values of c^c/b computed. If the discontinuity occurs at 
some other position, the value of (c^c/b)* must be interpolated. This 
interpolation may be accomplished in the following manner. Even though 
the spanwise lift distribution is continuous through this point, the 
point is singular. Its singularity is due to the singularity of the 
corresponding point of the distribution of c^c/bB which has an infinite 

slope but zero radius of curvature. At this point, from equation (11), 


C *2 C f rt0*sin 0* 
b8 } ~ 8l00 


( 12 ) 


Since both the distribution of cjc/b and the distribution of ci^c^bb con- 
tain the same type of singularity, a curve with c^c/b as the dependent 


variable and 


c 7 c 

l 2 


i 


bB as the independent variable would have no singular 
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point and could be approximated by a polynomial for which Lagrange's 
interpolation formula (reference 3) is applicable. From Lagrange's 
formula ■ 



(13) 


where 



(14) 


The number (n + 1) of terms retained determines the degree (n) of the 
polynomial used in the approximation. It has been found that using 
four terms (two on each side .of the desired point), which define a 
third-degree curve, gives values at the end of the flap which agree 
within about 1 percent with the values obtained by means of the method 
of reference 4 but with about 45 terms of a trigonometric series. 

Values of the interpolation multipliers determined in this manner are 
presented in table III for various spanwise positions of the end of a 
flap. This type of interpolation is suitable only if the variables are 
single -valued in the range of- the interpolation. This limitation pre- 
cludes the determination of multipliers for = q .1 but this position 

is out of the range of practical flap spans. 

The method of determining the lift distribution by means of succes- 
sive approximations is the same as that of reference 1 with the added 
step of determining the correction factor a c . For a given geometric 
angle of attack, the lift distribution is assumed. The uncorrected 
induced angle of attack is determined by equation (9) for asymmetrical 
distributions or by a corresponding equation using for symmetrical 

distributions. The multipliers 7 ^ for anti symmetrical distributions 

can be used only if the lift curves are linear and the value of 6 is 
independent of lift coefficient. • The value of (c z c/b)* is obtained by 
means of equation^(l3) and is divided by (c/b)* to determine c j*. At 
this value of c^ , 8 is read from the section lift curves as the. 
difference between the effective angle of attack for the section without 
the flap and that for the section with the flap, both at the station at 
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the end of the flap. This value of 5 is multiplied by the appropriate 
values of '<x c /& from table II to obtain clq . The values of a^. are 
added to the uncorrected values oty to obtain the values of induced 
angle of attack which are then subtracted from the geometric angle of 
attack to give the effective angle of attack at each spanwise station. 

At each value of effective angle of attack, the corresponding value 
of Cj is read from the appropriate section lift curve and multiplied 
by the value of c/b at that station to obtain a check value of cjc/b. 
If the check values do not agree with those originally assumed, a sec- 
ond approximation is made and the process is repeated. Further approxi- 
mations are made until one is found which is in agreement with the 
check values. 


Determination of Three-Dimensional Section Data 

I 

If the two-dimensional lift data were used directly, a discontinu- 
ity in the spanwise distribution of maximum lift coefficient^ would he 
indicated at the end of the deflected flap or aileron. Inasmuch as this 
discontinuity does not exist in three-dimensional flow, some means of 
altering the two-dimensional data must he used to obtain what may he 
called three-dimensional section data. One method of alteration is 
described in reference 5 for use with linear section lift data. In 
order to use this method, the- lift distribution must he broken up into 
the basic lift distribution due to twist and the additional lift dis- 
tribution due to angle of attack. When nonlinear section lift data are 
used, the lift distribution cannot be broken up in this manner and the 
method of reference 5 is not applicable. For this reason, the method 
hereinafter described has been developed. Although this development is 
not rigorous, the reasoning behind it is substantiated by experimental 
pressure distributions. 

The foundation of the method used herein is the calculation of two 
curves for which linear section lift data are used. One curve, the 
spanwise distribution of section lift coefficient for a wing with flaps 
but no other twist, is calculated, by the method described herein but 
adapted for linear lift curves in a manner similar to that 4 of refer- 
ence 1. The other curve is the additional lift-coefficient distribution 
(constant absolute angle -of -attack distribution) calculated as in refer- 
ence 1. These curves may be calculated for any convenient .angle of 
attack and for -any convenient value for the discontinuity in angle of 
attack at the end of the flap. A typical set of curves is shown in 
figure 5- It is reasoned that the root section of a wing with partial - 
span flaps would be acting most nearly like that* of* a wing with full- 
span flaps so that the additional lift-coefficient distribution is mul- 
tiplied by a suitable constant to give the same value at the root. 
Likewise, it is reasoned that the tip section would be acting most 
'nearly like that of a wing without flaps so that the additional 
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lift-coefficient distribution is multiplied by another constant to give 
the same value at the outermost station used in the computations (the 
0 . 9877 -semispan station when 10 points on the semispan are used). The 
differences between the values on the curves axe then divided by the 
difference between the values on the additional lift -coefficient curves 
at the end of the flap to obtain the factor F shown in figure 6 for 
several wings. In order to interpolate for the points on the additional 
lift-coefficient curves at the end of the flap, the multipliers p m , 
presented in table IV, have been determined. These multipliers were 
obtained by using 6 as the independent variable in equation (lk) 
instead of c^gC/bb. It can be readily seen that. the factor F is 

independent of the angle of attack and magnitude of the discontinuity 
used in the original computations. The factor F is also relatively 
independent of lift-curve slope so that it can be used near maximum lift 
where the lift curves axe nonlinear. 

The maximum lift- coefficient values axe altered by means of the 
factor F according to the relation 



The values of Cy and a axe then altered according to the equations 

^ altered 

c 



The edge-velocity factor E is used in the same manner as in refer- 
ence 1. The value of E given in reference 6 is, however, probably 
more, accurate than the ratio of semiperimeter to span used in refer- 
ence. 1 . From reference 6 for unswept wings, 


TcJJ 


unaltered 


‘•max 


fo 


‘'max 


(16) 



( 18 ) 


The foregoing description of the alteration to the two-dimensional 
section data has been limited to wings with symmetrical partial-span 
flaps. For wings with deflected ailerons, a similar method based upon 
the same reasoning could readily be devised. 
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The reasoning behind the method of altering the two-dimensional 
data is substantiated by the results shown in figures 7 and 8 and 
obtained from experimental pressure distributions for a wing of aspect 
ratio of 6.0, taper ratio of 0 .^>, and NACA 64-210 airfoil sections. 
Chordwise pressure distributions were obtained at. six spanwise stations. 
Two spans of flaps, 0.49b and 0.51b, were tested so that the pressure 

distribution at ^ = 0.5 was obtained, in one case, just outboard of 

the end of the flap and, in the other case, just inboard of the end of 
the flap. As shown in figure 7 , the pressure distributions at each sta- 
tion are very similar for the two flap spans, even at the station 0.5b/2. 
Furthermore, there is a gradual, spanwise change from the type of loading 
associated with an airfoil section with flaps to the type associated 
with a plain airfoil. 

For comparative purposes, the section normal-force coefficients, 
obtained by integration of the chordwise pressure distributions at each 
spanwise station, together with calculated values of section lift coef- 
ficient, interpolated for the same spanwise station, are shown in 
figure 8. The calculations for the flaps-neutral and full -span-flaps 
configuration were made according to reference 1 and those for the 
partial-span-flaps configuration were made as described herein. For 
the experimental data, jet-boundary and stream-angle corrections have 
been applied to the angle of attack, but no corrections have been applied 
to the values of c n for the effects of the model supports and the 
boom containing the pressure tubes. Although some of the disagreement 
shown may be attributed to the assumptions involved in the calculations, 
most of it is believed to be due to experimental inaccuracies both in 
the two-dimensional data used in the calculations and in -the three- 
dimensional data shown herein, inasmuch as similar disagreement is 
evident for the full-span configuration as for. the partial -span configu- 
ration. In general, however, the trends indicated in figure 8 serve to 
substantiate the method used to alter two-dimensional data. 

Wing Characteristics 

Numerical integration .- In the numerical integration described in 

reference 1, the multipliers n , n . 0 . and 0 were determined 

'm 7 'ms 7 m 7 ma 

by harmonic analysis. For example, it was shown that 



where 7^7 


sin 


— was designated as By inspection. 


it may be seen 
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that equation ( 19 ) could also be obtained through the application of 
the trapezoidal rule (reference 3) where f(2y/b) m sin (mit/r) is the 
function of uut/r and the interval between points is it/r . Although 
equation (l9).may be expected to give good results when f(2y/b) is 
approximately elliptic, in general the application of Simpson's 
(parabolic) rule to this problem has been found to give better results 
for a wider variety of curve shapes. Since r was originally assumed 
to be even, the interval m = 0 to r may be divided into r/2 regions, 
the area of each of which is 




sin 


(i + 1)* 
r 



sin 


(i + 2)« . 

r . 


(20) 


For the entire interval, m = 0 to r. 


= 2 g>(?) m v . (21) 

where tj is defined and used herein as 



The relations between ti . r\ , a , and a are the same as in 

nr ms 7 nr ma 

reference lj namely, N 


tj = 2t] (m f 

i) 

(23a) 

tj = tj fm = 

'ms 'm \ 

i) 

(23b) 

_ nut 

o m = — COS -r— 

m 2 r 


( 24 ) 
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a ma 2o m 


nut 

= ti cos — 
'm r 


( 25 ) 


Values of ti , ti , a , and a are given in table V for r = 20. 
m tus in mo. 

For cases for which the end of the flap or aileron is at a point 
where m is not an even integer, a special form of equation (20) is 
required which gives additional multipliers for numerical integration 
in the vicinity of the end of the flap or aileron. Such a case is illus- 
trated in. figure 9 and the individual areas are 


Area (1) = 


Area (2) = 


s(2s + 3t) „ 

6(1 VtT 1 


t3 ' 

6s(s + t) 1 


s(s + 3t) 


(26) 

6t f 2 

6(s + t)t f 3j 

t(t + 3s) r 
~Zs f 2 

. t(2t + 3s) -1 
+ S(s +~t) f 3j 

(27) 


Area (l) + (2) = — 


(2s - t)(s + t) 

6s 1 


(s + t )3 f (2t - 
+ &t 2 


s)(s + t) 

Zt f 3 


(28) 


where s and t are expressed as fractions of Jt/r. It can be seen 
that' equation (28) is the same as expression (20) when s = t =-l. An 
example of the computation of these special multipliers is given in 
table VI for 0.5-span flaps. For the area B in the sketch accompanying 


this table, equation (28) is used with 
area C, equation (26) is used with s = 


s = 1 and t = 
~ and t = 1. 


and for the 
Values of the 


multipliers are given in table VII for various locations of the end 
of a flap. Values of the multipliers q , cx m , and can be readily 
obtained from the values of t| ms by means of equations (23) to (25). 


It should be noted that two multipliers are given at the end of the 
flap. For distributions which are continuous through the end of the 
flap, such as the distribution of cjc/b, the sum of the two multipliers 
may be used. For distributions which are discontinuous at this point, 
such as the distribution of c^^^/b, each multiplier should be used 

separately with its appropriate value. For example, for = 0.5, the 
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multiplier -0.03023 should be used with the value at the outboard end 
of the flap, and the multiplier 0.03930 should be used with the value 
at the inboard end of the unflapped span. 

Wing coefficients .- The formulas derived in reference 1 are 
repeated herein for convenience. For asymmetrical distributions 



For symmetrical distributions 



For asymmetrical distributions 



For symmetrical distributions 

p JtA £^ / c Z c \ 


(29a) 


(29b) 


(30a) 


(30b) 


For asymmetrical distributions 



For anti symmetrical distributions 



For asymmetrical distributions 



(31a) 


(31b) 


( 32 ) 
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The values of ti , n , a . and a should be used as defined herein, 
'm' ins' nr ma 7 

and the summations should also include the special, multipliers in the 

vicinity of the end of the flap or aileron. 


ILLUSTRATIVE EXAMPLES 
Symmetrical Distributions 


The method described is applied herein to a wing, the geometric 
characteristics of which are given in table VIII. The section lift 
curves are shown in figure 10 as dotted lines before being altered and 
as solid lines after being altered. The unaltered data are obtained by 
interpolation of two-dimensional data for the proper Reynolds number and 
airfoil thickness -chord ratio for each spanvise station in a manner 
similar to that of reference 1. The altered curves are obtained as 
described previously and are used in the calculation of the lift dis- 
tributions in tables IX and X. These tables were designed for use with 
calculating machines Capable of performing accumulative multiplication. 
The; mechanics of computing are explained -in the tables, but a few items 
need additional explanation. 

The initial approximation of c^c/b is obtained in table IX. The 
computations in this table are based on the method of reference 6, the 
pertinent equations of which are modified to suit the present purpose. 
The spanwise stations, at which values of c^c/b are computed, are 
listed in column (l). Columns (2) to (4) are used to obtain the addi- 
tional. lift distribution c jc for Cl = 1 according to the 

approximate relation 



Columns (5) to (7) are used to obtain faired values of the twist €5 
due to flap deflection according to the relations 
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(34b) 


This fairing is illustrated in figure 11. Columns (8) and (9) are used 
to obtain the average twist according to the equation 



(35) 


The value of obtained by this summation is given at the bottom of 
the table. Columns (10) to (12) are used to obtain an angle of attack 
a' defined, for the flapped span, as 


AE 

AE 


7“§( e 6 " € s) + *6 _ 5 


+ a 


(36a) 


and, for the unflapped span, as 

a ’ = inri( c & - r &) + + a (36b > 

The values of the geometric angle of attack a in column (11 ) include 
the values of continuous twist €. Values of c^ a t^ in column (13) 

are obtained for the section lift curves at the angle a 1 . In t£e non- 
linear range of the section lift curves, different values of Cj ( a * ) 

are found corresponding to the two values of a' at the end of the 
flaps. For the purpose of this table, the average value of c ^*( a ! ) 

used. Finally, columns (lU) and (15) are used to obtain the initial 
approximate values of c^c/b according to the relation 



approx 


c 7 c 
z al 

c(AE + 2) c *(a’ ) 


which was adapted from the equation 


C Z C _ a «o C *al C 

_ AE 


c\ 

b ~ AE + 2 c 

e + AE 

+ £ \ 

- V 


(37) 


( 38 ) 


given in reference 6. 
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T^e value of & used in table IX should correspond to the value 
of Cj obtained at the bottom of the table. Some value of 6 must 
be selected, ^however , for the computations needed to determine the 
value of C 7 . It is therefore advisable to omit the computations in 
columns (10) to ( 15 ), except the bottom row, until a value of 6 has 
been found by trial and error which is consistent with the value of c \ . 

If the lift distribution is to be obtained at a relatively high 
angle of attack, as in the example shown, the values of a' may be in 
the rapidly curving part of the section lift curves or even greater than 
the angle of attack for maximum section lift coefficient-. For the pur- 
pose of this table, the section lift curves may be extended in the gen- 
eral direction of the nearly linear part of the curves in order to 
obtain the values of c^( a ')* This procedure is justified inasmuch as 

the subsequent operation in column ( 15 ) will reduce the value of cj. 

The approximate values of c^c/b (column (15) , table DC) are used 
in column (3) of table X and the computations indicated are performed. 
Except for the computation pf a c , the computations in this table are 
the same as those in reference 1 for wings without flaps. The value 
of 5 used in column (l6) is obtained at the section lift coefficient 
Cj corresponding to the value of (c^c/b)* obtained by means of the 
interpolation multipliers v m . 

In this method of successive approximations, the value of c^c/b 
in column (22) will usually not check the initial approximate values in 
column (3)* The values to be used in subsequent approximations may be 
found by the equation 



where A'(c^c/b) is the increment to be added to the approximate values 
to obtain- succeeding approximate values, A(c£c/b) is the difference 
between the check values and the approximate values (column (22) minus 
column (3)), and K and Ki are constants for any particular wing. 
Equation (39) is derived in the appendix, and values of K and Kj. 
for r = 20 are presented in figure 12 as functions of AE/q. These 
values compare favorably with those empirically determined in refer- 
ence 1 (Kp - 3 , % = 1, = 0 , and K = 8 to 10 j . Although these 

values were obtained for elliptic wings, they can be used for wings of 
other plan form. The number of terms of equation (39) -needed for any 
particular approximation depends upon the convergence of the approxima- 
tion; fewer terms are needed when the differences A(c 2 c/b) are smalls or 
when positive differences nearly cancel negative differences. Values 
of. Ki for values of i greater than 3 are small enough to be considered 
negligible . ' 
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Some additional explanation of the use of equation (39) may be 
desirable. For values of m from 3 to 7 , equation (39) may be used 
directly if the constants ^>3 are neglected. For values of- m 

from 8 to 10, equation (39) may be expanded as 



(toe) 


for symmetrical distributions, since a(^J = . For values 

' • J ' 'r-j 

of m of 1 and 2, equation (39) may be expanded as 
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since 



After convergence is obtained in table X, the values of c^c/b 
(column (3)) and ccq (column (l8)) are entered in table XI and the lift 
and induced-drag coefficients axe determined through use of the appro- 
priate multipliers t> . • 


Asymmetrical Distributions 

If the angle -of -attack distribution is not symmetrical, the asym- 
metrical multipliers must be used in the nonlinear range of the 

section lift curves. Typical asymmetrical distributions are those for 
a rolling wing or for a wing with deflected ailerons. Illustrated in 
table XII is the case of a wing with flaps but without deflected ailerons, 
which is rolling at such a rate that the tip helix angle pb/2V gener- 
ated is 0.01 radian or 0.573°- Added to the normal angle of attack 

is an increment equal to ^ (in deg), which is the equivalent 

twist of the rolling wing. In order to reduce the size of the computing 
form, the table of the multipliers (3^ is arranged in the form shown. 

The values of (3 m ^ for a positive value of 2y/b are the reverse of 
those for the corresponding negative value. Instead of reversing the 
values of 3^ in table XII, the values of cjc/b are written in 
reverse order in column (14). Using these values with the values of 
gives the uncorrected angle 0^ for the stations listed at the 
bottom of the table. ■ . 

In general, the value of 6 will be different for the two sides 
of the wing, so that the appropriate values must be used with columns (16) 

and (17 ) • The values of a c /S_ in column (l6) are those for = _o.6 

from table II, while the values of a c /& + in column (17) are the 

2v 

negative of those for =0.6. 

Another modification must be made for asymmetrical lift distribu- 
tions since the edge-velocity factor E' should be used for the antisym- 
metrical part of the distribution (see reference 6), whereas the edge- 
velocity factor E has been used to alter the two-dimensional lift 
curves. From reference 6 for unswept wings, 

E ’ -\F+M ■ 
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This value may be taken into account in the following manner: The 

symmetrical part of the effective angle distribution is 

(a - ai) k + (a - <n) p _ k 
- _ 

which is used directly with the altered lift curves. The antisymmetrical 
increment in the angle distribution is 


(a - a ± ) - (a - 04) 

± k x r-k 

which must be multiplied by the ratio E/E' in order for it to be used 
with the same lift' curve. The effective angle is therefore 

(a - 04 ) k + (a - C4) r _ k E (a - a i ) k - (a - 04)^ 

“ 2 + E 7 2 : 


= (»• - »l) k - ^ - «!>* - (a - Oi) r J 

= a - a ± - (43) 

where 

= E '2E' E [^ a " " (a “ a i) r .k] (U) 

Equation (44) is confuted in column (21) of table XII. 

Other than these modifications, the computing required for table XII 
is similar to that for table X. 


DISCUSSION 


The lift characteristics of two wings without flaps and with 
60-percent and full-span flaps have been calculated by the method 
described herein and are presented in figure 13 together with experi- 
mental results from reference 7. One wing had NACA 64-210 sections and 
was equipped with split flaps. The other wing had NACA 65-210 sections 
and was equipped with split, single slotted, and double slotted flaps. 

For the split-flap conditions, the agreement between calculated an d 
experimental results is quite satisfactory, whereas the agreement is less 



NACA TN 2283 


23 


satisfactory for the single- and double -slotted-f lap conditions. Since 
the discrepancies occur for both the 60-percent-span and full-span con- 
ditions, they are probably due to differences between the two-dimensional 
and three-dimensional section characteristics rather than due to the 
method of calculating. Some of the discrepancies in maximum lift coef- 
ficient may be due to the fact that the characteristics of these wings 
were extremely sensitive to small surface irregularities. 

The stalling characteristics of these same wing-flap combinations 
are presented in figure Ik together with the calculated stall-margin 
distributions. The stall margin Ac j is the difference between the 
maximum section lift coefficient altered as described herein and the 
section lift coefficient at the maximum wing lift coefficient. The 
spanwise location of zero margin should correspond to the location of 
initial stall, and the margin at other spanwise locations is an indica- 
tion of the manner in which the stall spreads. In general, the agree- 
ment between the experimental and calculated stalling characteristics 
is very good. 

The foregoing comparisons between calculated and experimental 
results were made for the same Reynolds number. If possible, such com- 
parisons should be made at the same Mach number also, unless the Mach 
number is low enough to have a negligible effect. Even at relatively 
low values of free-stream Mach number, adverse compressibility effects 
on maximum lift coefficient have been noted (reference 8) when sonic 
velocity is reached locally on a wing. Similar effects in two- 
dimensional flow have not as yet been thoroughly investigated so that 
calculations based on available two-dimensional data must.be limited to 
subcritical Mach numbers. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va. , November 13, 1950 
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APPENDIX 


DETERMINATION OF COEFFICIENTS' USED TO 
OBTAIN SUCCEEDING APPROXIMATIONS 


In the method of successive approximations to determine the lift 
distribution, it is desirable to reach convergence with a minimum number 
of approximations. This desideration necessitates that each successive 
approximation be obtained from preceding computations . in some maimer. 

The manner in which these operations were performed in reference 1 and 
the coefficients used therein were determined empirically. It is shown 
hereinafter, however, that the procedure and similar coefficients can 
be determined theoretically. 

In the following derivation, A(c^c/b) is used to designate the 
difference between the check values and the approximate values, and 
A'(cjc/b) is used to designate the increment to be added to the approxi- 
mate values to obtain the succeeding approximate values. The section 
lift curves are assumed to be linear and 6 to be zero. The operations 
performed during the first approximation (table; X) may then be repre- 
sented by the equation , 



If suitable increments A'(c^c/b) are chosen so that the check values 
become equal to the approximate values for the second approximation. 



The difference between these equations is 




(A4) 


From reference 1 
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The lift -curve slope may he expressed as 

.2 


T] Jt 

a o E 90 


(A5) 


where r\ is the ratio of the actual two-dimensional lift -curve slope 
to the theoretical thin-airfoil value. For elliptic wings y 


kit 


'c\ _ 4 , 

(b)k JtA S n r 


(a 6) 


Therefore 


a o c \ _he. 

b ). - AE 


(AT) 


which, for constant r \ , is constant for all spanwise stations, 
stituting this value into equation (A3) yields 


Sub- 




+ 

k m^k 


^mk 

^-kk 




Equation (A8) represents r/2 simultaneous equations which may be 
represented in matrix form as 



r , c z c_ 1 

ae 

c z c 

M 

A ' — 

■ qr 

L b J 


(A9) 


where [jdJ is a matrix with all the principal diagonal elements equal 

to 1 + — and the other elements are \nkAkk* This matrix can 

readily be put into a symmetrical form and its- reciprocal obtained by 
one of the standard methods presented in reference 9- Then 



(A10) 


For a given value of AE/t], equation (A10) may be expressed in the form 



(All) 
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because of the particular properties of this reciprocal matrix. For 
convenience Kp is made equal to unity and the values of K and Kq 
are adjusted accordingly. The values given in table XIII and figure 12 
were obtained for various values of AE/t] and r = 20. . 
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TABLE V.- WING-COEFFICIENT MULTIPLIERS 
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TABLE VI.- CALCULATION OP AREA MULTIPLIERS WITH SPECIAL 
CONSIDERATION AT THE END OP A 0 . 5 ^SPAN FLAP 


ii 


b 

m 

0 

10 

.156 

9 

.309 

e 

.454 

H 

.500 

6 I 

CO 

00 

LTN 

• 

6 

.707 

5 


Factors for areas indicated 
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TABLE VII.- AREA MULTIPLIERS FOR VARIOUS SPANWISE LOCATIONS OF THE END OF A FLAP 
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TABLE DC.- CALCULATION OF INITIAL APPROXIMATION OF LIFT DISTRIBUTION FOR EXAMPLE WING WITH 0.6 SPAN FLAPS 































TABLE X — CALCULATION OF LIFT DISTRIBUTION FOR EXAMPLE _ WING WITH 06 SPAN FLAPS 
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TABLE XE— CALCULATION OF ASYMMETRICAL LIFT DISTRIBUTION FOR EXAMPLE WING WITH — CLfi SPAN FLAPS 
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TABLE XIII.- COEFFICIENTS USED TO OBTAIN 
SUCCEEDING APPROXIMATIONS 


AE 

n 

K 

Ko 

K 1 

K2 

k 3 

4 

10.095 

2.384 

1.000 

0.498 

0.339 

8 

8.418 

3.000 

1.000 

, 401 

.258 

12 

' 8.141 

3.571. 

1.000 

• 338 

.215 
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Figure 1.- Typical lift distributions for a wing with 0.5-span flaps. 









Figure 4.- Values of the lift distribution per degree of discontinuity in 
induced angle of attack as a function of the spanwise location of the 
discontinuity 2y*/b. 
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Figure 5.- Typical curves used to obtain factor for altering two-dimensional 
. . data. 



Figure 6.- Factor for altering two-dimensional data for several wings. 
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(c) NACA 65-210 sections; single slotted flaps. 


Figure 13*- Continued. 
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(d) NACA 65-210 sections; double slotted flaps. 
Figure 13.- Concluded. 
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(d) NACA 65-210 sections; double slotted flaps 








